AD- A  1 20  996 


A  UNIQUENESS  THEOREM  IN  ELASTODYNAMICS<U)  WISCONSIN 
UNI V-MAD I  SON  MATHEMATICS  RESEARCH  CENTER  D  S  JONES 
SEP  82  MRCTSR-2425  DAAG29 -80-C -004 1 


UNCLASSIFIED 


F/G  20/ B 


l 

I 

I 

I 


i 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-  1963-A 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison.  Wisconsin  53706 

September  1982 

(Received  August  23,  1982) 


DTIC  FILE  COPY 


Appravad  far  paklic  ralaasa 
Distrikatiaa  anliaritad 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


8*  11  02  06  4 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 


A  UNIQUENESS  THEOREM  IN  EIASTODYNAMICS 
D.  S.  Jones 

Technical  Summary  Report  #2425 
September  1982 


ABSTRACT 


A  uniqueness  theorem  is 
elastic  waves  by  a  body  with 


established  for  the  scattering  of  harmonic 
continuously  varying  parameters  placed  in  a 


homogeneous  medium. 
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□  □ 


SIGNIFICANCE  AND  EXPLANATION 


The  passage  of  waves  through  a  solid  is  a  natter  of  sons  practical 
importance  because  of  their  effect  on  structures  and  also  because  they  can  be 
responsible  for  the  transmission  of  noise.  Vibration  may  occur  naturally  due 
to  such  thing  as  earth  movement  and  wind  gusting  or  may  be  caused  by 
machinery.  Sometimes  oscillations  are  deliberately  induced,  as  in  ultrasonic 
testing,  to  check  the  strength  of  bonds,  to  detect  flaws  or  to  locate  pockets 
of  material  different  from  their  surroundings.  Theoretical  prediction  of 
phenomena  is  therefore  a  vital  adjunct  of  investigation  in  the  field. 

The  theory  is  baaed  on  computing  solutions  of  the  equations  of  a 
mathematical  model.  These  equations  may,  and  usually  do,  have  several 
solutions  so  there  is  a  task  of  identification  and  interpretation  of  their 
relevance)  clearly,  wasted  effort  will  be  avoided  if  it  can  be  indicated  which 
solution  or  solutions  to  seek.  Previously,  it  has  been  shown  that,  if  a  loss- 
free  object  of  constant  properties  is  embedded  in  another  substance  with 
constant  material  properties,  there  is  only  one  solution  which  allows  energy 
to  radiate  away  from  the  body  after  it  has  been  struck  by  an  incoming  wave. 
This  report  extends  the  theory  to  obstacles  whose  properties  change 
continuously  from  point  to  point  and  also  permits  them  to  have  losses.  It 
shows  that  there  is  still  a  single  solution  to  be  found  provided  that  one 
keeps  to  the  mathematical  criterion  that  the  disturbance  from  the  obstacle 


spreads  outwards 


A  UNIQUENESS  THEOREM  IN  ELASTODVNAMICS 
D«  S.  Jones 

1.  Introduction 

In  his  book  on  elastodynamics  Hudson  (1980)  points  out  the  need  for  a 
uniqueness  theorem  for  materials  which  are  not  homogeneous.  This  paper  makes 
a  contribution  to  filling  that  gap  by  providing  a  uniqueness  theorem  for 
inhomogeneous,  but  isotropic,  bodies.  Inhomogeneities  of  two  types  will  be 
examined.  Those  in  which  the  material  parameters  are  constant  except  across 
certain  surfaces  of  discontinuities  and  those  in  which  the  parameters  vary 
from  point  to  point  but  have  some  continuity  and  differentiability  available. 
The  piecewise  homogeneous  case  has  already  been  discussed  by  Kupradze  (1963) 
but  it  is  included  here  partly  for  completeness  and  partly  because  it  involves 
only  a  short  argument  from  formulae  which  are  needed  for  the  continuous 
inhomogeneity.  In  our  investigation,  the  material  parameters  are  permitted  to 
be  complex,  with  some  limitations,  so  that  lossy  substances  are  not  excluded 
from  the  theory,  only  bonded  bodies  will  be  the  subject  of  study. 

Section  2  sets  out  the  basic  linear  equations  which  are  of  concern  and 
formulates  some  of  the  constraints  on  the  field,  standard  formulae  for  the 
representation  of  the  displacement  in  a  homogeneous  body,  whether  finite  or 
infinite,  are  given  in  Section  3. 

As  it  turns  out  a  proof  of  uniqueness  for  the  general  problem  revolves 
about  a  proof  for  the  interior  of  a  finite  body.  This  interior  question  is 
examined  for  the  piecewise  homogeneous  substance  in  Section  4  and  for  the 
continuous  inhomogeneity  in  Saction  5.  While  the  analysis  of  Section  4  is  a 
straightforward  application  of  the  representation  for  the  displacement  that  of 
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Section  5  is  much  more  complicated  and  probably  the  most  difficult  part  of  the 
whole  exercise.  Nevertheless,  the  theory  here  has  no  restrictions  (other  than 
continuity  and  differentiability)  on  the  complex  values  of  the  aaterial 
parameters  though  this  freedom  has  to  be  abandoned  in  the  full  problem 
later*  By  Mans  of  the  interior  theorems  the  relevant  theorems  for  the 
infinite  medium  are  derived  in  section  6  and  their  application  to  the 
uniqueness  problem  in  scattering  indicated  in  Section  7. 

Two  appendices  contain  results  whose  derivation  would  have  interrupted 
the  flow  of  argument  in  the  main  text.  Appendix  A  gives,  for  easy  reference, 
certain  properties  of  spherical  harmonics  needed  in  Section  5.  Appendix  B 
covers  the  basic  uniqueness  theory  and  expansion  properties  in  an  infinite 
homogeneous  medium,  in  particular,  it  verifies  that  one  of  the  two  customary 
radiation  conditions  can  be  disposed  of  without  affecting  uniqueness. 

Equations  from  the  Appendices  when  referred  to  in  the  main  text  are 
distinguished  by  the  appropriate  letter. 

t 

2.  The  governing  equations 

The  first  problem  to  be  considered  is  that  of  harmonic  elastic  waves  in 
an  isotropic  body  of  finite  sise.  The  body  occupies  the  voIum  T_  and  its 
surface  will  be  denoted  by  s.  The  voIum  outside  S  will  be  identified  by 
T+.  In  T_  the  displaceMnt  at  the  point  x  is  u(x)  with  Cartesian 
components  Uj  and  the  stress  tensor  is  T^.  The  Mtsrial  occupying  T_  is 
of  density  p  and  its  elastic  properties  are  specified  by  the  Lame  parameters 
X,  u.  It  will  be  assumed  that  there  are  no  body  forces.  Then  the  equations 
to  be  satisfied  in  T_  are,  when  the  tiM  dependence  is  e  , 


■2- 


(2) 


2 

u  pu. 
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where  6 ..  is  the  standard  Kronecker  symbol  and  the  usual  summation 
jk 

convention  has  been  enployed.  me  Cartesian  for*  of  the  governing  equations 
is  given  with  (x1,x2,x3)  the  coordinates  of  the  point  x. 

In  eost  applications  p  is  positive  and,  in  lossfree  media,  X  >  0  and 
W  >  0.  When  dissipation  is  present  X  and  u  can  have  imaginary  parts.  For 
much  of  the  subsequent  analysis  they  can  be  taken  as  arbitrary  complex 
quantities  but  it  will  always  be  assumed  that,  on  T_  and  S, 

I Vl  >  1/K#  I X+2ul  >  1/K  »  (3) 

where  K  is  a  positive  finite  constant. 

The  solutions  of  (1)  and  (2)  to  be  found  depend  on  the  boundary 
conditions  on  s  and  whether  there  is  any  transfer  of  energy  between  T_ 
and  T+.  When  wave  motion  in  T+  has  to  be  taken  into  account  the  media  will 
be  assumed  to  be  bonded  across  the  interface  s  so  that  the  traction  and 
displacement  are  continuous  there. 


3.  The  homogeneous  medium 

For  subsequent  purposes  it  will  be  convenient  to  have  a  representation  of 
the  displacement  when  the  medium  is  homogeneous  with  X,  |i,  p  having  the 
constant  values  \Q,  uQ,  pQ  respectively.  It  will  be  supposed  that  ?  0 
in  concordance  with  (3).  Let 

,  .  1  a2  e-i<*lgrxl 

gjk<s*x  4wPoW2  41%isxi  V 

where  a  ■  {Pq/(Xo+2Uq))^  *  b  ■  IPq/^q)^2  &ra  the  slownesses  of  F-  and  8- 
waves  respectively.  For  complex  parameters  the  square  roots  are  defined  by 
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their  principal  values.  Zn  view  of  (3)  both  a  and  b  are  finite.  The 
ranges  of  the  parameters  will  be  restricted  so  that 


Ia(ua)  <  0,  Xa(cifc)  <  0  (5) 

but  the  real  and  imaginary  parts  of  ua  are  not  permitted  to  vanish 
simultaneously  nor  are  those  of  ub.  If  Im(ua)  -  0  then  ua  must  be 
positive  and  similarly  for  ub. 

The  tensor  g^  specified  by  (4)  satisfies 


a2g  32g 

( W  15^  +  M0  9^  *  po“2gik  "  5ik6(5“X> 


(6) 


By  means  of  the  divergence  theorem  and  (1),  (2),  (6),  Betti's  representation 
(see,  for  example,  Kupradse  (1963))  can  be  obtained,  namely 


V*>  ’  VVa  V*'*’  *  V "iV-jV  ^  Vs'*’ 


(7) 


for  x  in  T_.  In  (7)  n  is  the  unit  normal  to  8  out  of  T_  into  T+. 

If  x  is  in  T+  the  left-hand  side  of  (7)  is  replaced  by  zero. 

Suppose  now  that  T+  is  a  homogeneous  isotropic  medium.  A  representa¬ 
tion  similar  to  (7)  can  be  derived  provided  that  appropriate  radiation 
characteristics  at  infinity  are  prescribed.  There  are  two  possibilities  that 
will  be  considered. 

In  the  first  case  suppose  that  Im(ua)  <0  and  Im(ufc)  <  0.  It  is  then 
evident  from  (4)  that  g^y.  decays  exponentially  at  infinity.  Assume  that  the 
displacement  and  stress  also  are  exponentially  attenuated  at  infinity.  Then 
appiy  (7),  which  is  valid  for  the  interior  of  any  closed  surface,  to  the 


volume  between  a  large  spherical  surface  and  the  exterior  side  of  S.  Because 


of  the  exponential  falling  off  of  g^,  u,  the  contribution  of  the 

spherical  surface  tends  to  zero  as  the  radius  of  the  sphere  approaches 

infinity*  Therefore,  if  the  definition  of  n  is  retained  unaltered,  u^(x) 

is  the  negative  of  the  surface  integral  in  (7)  when  x  is  in  T+.  Of  course, 

it  must  be  remembered  that  in  this  surface  integral  u  and  T,.  have  the 

~  JX 

values  on  S  which  are  approached  from  T+  whereas  when  the  integral 
represents  u^  in  T_  the  values  of  u  and  on  the  interior  side  of 

s  must  be  employed. 

The  second  case  is  that  in  which  lm(ua)  ”  0  and  X®( ub)  ■  0*  Bow  it  is 
assumed  that,  as  R  m  |x|  ♦  *»,  Ru  is  bounded  and  that 

R{Vjm  +  iUbW0(ttj  '  ".Vj1  +  iwa(  V2VXjUmXm}  *  0  (8) 

A 

where  x  is  a  unit  vector  in  the  direction  of  x.  Actually,  it  is  shown  in 
Appendix  B  that  the  requirement  for  Ru  to  be  bounded  is  superfluous  but  it 
is  easier  to  justify  some  of  the  subsequent  statements  if  this  condition  is 
retained.  Bearing  in  mind  that,  as  |x|  ♦  •  with  £  fixed. 


g  (  ,  ~  [b2  -iU*>(R-X.£)  _  2  -iwa(Rrx^) 

9jk  4*pftR  1 


(9) 


4*U0R 


k_  e-iui>(R-x*£) 


we  see  that  the  contribution  from  the  sphere  at  infinity  again  vanishes.  The 
conclusion  is  that  can  be  represented  in  T+  by  the  negative  of  the 

surface  integral  in  (7). 


Zt  is  conceivable  that  one  of  Zm(ua)  and  Zm(tito)  could  be  zero  while 
the  other  was  negative  but  this  does  not  seem  to  be  of  sufficient  practical 
significance  to  justify  separate  consideration.  Therefore,  only  the  two  cases 
already  mentioned  will  be  discussed  subsequently. 

4.  The  interior  problem  for  piecewise  homogeneity 

This  section  and  its  successor  are  devoted  to  examining  what  can  be  said 
about  the  field  in  T_  when  the  displacement  and  traction  vanish  on  S. 
Obviously,  if  the  medium  is  homogeneous  (7)  may  be  invoked  with  the  conclusion 
that  u  =  0  in  T_.  Zt  follows  from  (1)  that  the  stress  also  disappears. 

Actually  (7)  can  be  extended  to  piecewise  homogeneous  bodies  in  suitable 
circumstances.  Let  S1  be  a  closed  surface,  entirely  inside  S,  where  the 
material  properties  change  discontinuous ly  from  one  set  of  constant  values  to 
another  set.  Apply  (7)  to  the  closed  surface  consisting  of  S  +  S1 .  The 
integral  over  S  is  removed  by  the  boundary  conditions.  As  for  the  integral 
over  it  represents  a  field  everywhere  outside  S1  and  it  is  obviously 

analytic.  But  it  is  identically  zero  outside  S  and  therefore  analytic 
continuation  ensures  that  it  is  identically  zero  outside  S1 .  Therefore  the 
traction  and  the  displacement  on  the  exterior  of  are  zero.  Zf  the  media 

are  bonded  across  the  traction  and  displacement  are  continuous  and 

therefore  zero  on  the  inner  side  of  .  Zf  the  interior  of  is 

/ 

homogeneous  the  preceding  paragraph  tells  us  that  the  field  is  zero.  Zf  there 
is  another  closed  surface  S2  inside  S.,  where  the  material  properties 
change  but  in  a  bonded  manner  we  repeat  the  argument  starting  from  the  closed 
surface  S1  +  S2<  Obviously  the  procedure  may  be  continued  for  any  finite 
number  of  closed  interfaces  each  totally  enclosed  by  its  predecessor.  Such  an 
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arrangement  nay  be  distinguished  by  the  adjective  nested*  There  is,  of 
course,  no  reason  why  Sj  should  not  consist  of  a  finite  number  of  distinct 
closed  parts. 

It  oust  be  emphasized  that  the  above  process  may  break  down  if,  at  any  of 
the  interfaces,  the  boundary  conditions  are  other  than  the  continuity  of 
traction  and  displacement.  In  order  to  stress  this  point  the  word  bonde*  is 
included  in  the  following  theorem  which  has  now  been  demonstrated. 

Theorem  1 .  if  the  traction  and  displacement  vanish  on  the  bounded 
surface  S  of  a  bonded  nested  piecewise  homogeneous  body  T_  the 
displacement  and  stress  tensor  are  identically  zero  in  T_. 


5.  The  interior  problem  for  continuously  inhomogeneous  bodies 

While  Theorem  1  does  cover  certain  inhomogeneous  bodies  the  departure 
from  homogeneity  consists  essentially  of  discontinuous  changes  across 
surfaces.  In  this  section  the  material  parameters  will  be  assumed  to  vary 
from  point  to  point  but  discontinuities  will  not  be  penitted.  In  fact,  it 
will  be  assumed  that  p  is  continuous,  \  is  continuously  differentiable  and 
that  |i  is  twice  continuously  differentiable. 

Throughout  solutions  of  (1)  and  (2)  in  T_  will  be  sought  in  which  u 
and  3u^/3x^  are  continuous.  Furthermore,  the  boundary  conditions 


u  ■  0,  n.T 

<W  #Wr  1  ' 


3  jk 


(10) 


will  be  imposed  on  S. 

If  T_  were  a  piecewise  homogeneous  medium  the  boundary  conditions  (10) 
would  be  sufficient  to  ensure  that  the  field  vanished  identically  but  the 
argument  of  the  preceding  section  cannot  be  carried  over  to  the  case  when  the 


7 


material  parameters  are  continuously  variable  because  the  representation  (7) 
is  no  longer  available.  Instead  one  is  forced  to  proceed  in  a  more  indirect 
fashion. 

In  order  to  fix  ideas  the  value  zero  is  assigned  to  1,11,0  in  T+. 


field  ul1). 

T0) 

jit 

(1) 

M) 

u  *  u. 

jk 

(1)  « 

.(1) 

u  *  0 , 

IW  »v 

T 

jk 

jk 


(x  e  T_ )  , 

<*/*  T_)  • 


Then,  under  the  assumed  conditions.  *  duj^/^x^,  x*j^  ,  3tj^^xk  are 

continuous  everywhere. 

Proof.  Until  the  stated  continuity  has  been  established  values  on  S  when 
x  approaches  there  from  T_  will  be  denoted  by  (  )_  and  those  for  x 
starting  in  T+  will  be  signified  by  (  )+. 

Continuity  in  T+  is  immediate  because  all  quantities  vanish  and, 
moreover, 

(u(1))+  "  Ou51)/»xJ{)+  -  0,  (T^>+  -  0,  0^)/3xk)+  -  0  .  (It) 

In  T_,  u(1)  is  continuous  by  assumption  and  (10)  gives  (u*1  * )_  »  0. 
Thus  the  continuity  of  u^1 ^  has  been  verified.  It  is  also  clear  from  the 
hypotheses,  (1)  and  (2)  that  Suj1  ^/3x^,  xjj^  and  are  continuous 

in  T_.  Accordingly,  the  lemma  will  be  confirmed  once  the  requisite 
continuity  across  S  has  been  demonstrated. 

Choose  any  closed  circuit  C  on  S.  Then,  by  Stokes's  theorem, 

/c(«[1)).  dg,  “  /s,  ha  (grad  u*1})_  dS  (12) 

where  S'  is  the  portion  of  S  within  C.  By  assumption  (10),  the  left-hand 
side  of  (12)  is  zero.  Since  C  is  arbitrary,  the  conclusion  is  that 
n  A  (grad  uj^)_  is  zero  on  S.  A  vector  product  with  n  then  enforces 
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on  S« 


Also  (10)  says  that  n.  *  0  on  S  and  so,  from  (1), 

j  JK  “ 


a  (1)  a  (D 
3u  3^ 


Xnk(-5iH-  +  +  T^- 


The  insertion  of  (13)  leads  to 


3um  3u"'  3u/ 1 ' 

av«(-3T~)-  +  •“W-siH- +  "( *H- 


on  S.  Multiply  (14)  by  n^  and  then,  by  virtue  of  (3), 

..j(3uJ1)/3n)_  -  0  .  (15) 

On  the  other  hand,  if  is  the  alternate  tensor,  multiplication  of  (14) 

bY  Apqknq  9ivea 

A~,vn„<3ui1>/3n>  *  0  (16) 

pqk  q  k 

on  invoking  (3)  and  noting  that  A^^n^n^  “0.  Equations  (15)  and  (16)  may  be 

expressed  as  n#3u^V3n  *  0,  nA  3u^V3n  ”  0  and  therefore  Ou^Vdn)  “  0 

on  S,  One  infers  from  (13)  that  { Su^1  ^/3x_j )_  *  0  on  S.  The  required 

continuity  of  3u^  /3x^  now  follows  from  (11). 

In  view  of  (1)  and  (11)  the  continuity  of  may  be  inferred  whereas 

that  of  a  con8e<luence  (2),  (11)  and  what  has  already  been 

established  about  u^1 ^ .  The  lemma  is  proved. 

Theorem  2.  In  T_  let  u,  3u^/3x^  be  continuous  and  satisfy  ( 1 ) ,  ( 2 ) 

under  the  given  assumptions.  If  u  ■  0  and  n .  t ..  »  0  on  s  then  u  5  0. 

1  -  ~  ~  -  j  jk  —  -  *»  ' 

T ..  =0  in  T  . 


are  first 


Proof.  To  permit  analysis  in  the  whole  of  space  u  and 
extended  to  u^1  *  and  as  in  Lemma  1.  However,  for  simplicity  of 

writing  the  affix  (1)  will  be  dropped,  on  the  understanding  that  its  presence 
is  allowed  for. 

Pick  as  origin  any  point  which  can  serve  as  the  centre  of  sone  sphere 
within  which  the  field  vanishes  identically.  Clearly,  there  are  many  possible 
choices  since  any  point  of  T+  has  the  desired  property. 

Let  R,  6,  $  be  spherical  polar  coordinates  based  on  this  origin.  There 
are  2n  +  1  independent  surface  harmonics  of  order  n.  Construct  an 
orthonormal  set  from  them  and  let  a  typical  member  of  that  set  be  Sn^(  6,  +) 
with  j  ■  -n,-n+1,...,n.  Then 

v2{Rnsnj(e,$)>  -  0  , 

v2{R_n“1sw.(e,^)}  -  o  <r  >  o)  , 

nj 

/_  S  J  .  dft  -  S .. 

■'$)  nj  nk  jk 

where  Q  is  the  surface  of  the  unit  sphere. 

Next,  define 

R  0 

On  account  of  (17)  and  (18) 

V2*nj(o,R)  -  0  (21) 

for  R  >  0.  in  addition  i|>  .(0,0)  *  0. 

nj 

The  function  ^  has  a  singularity  at  the  origin  but,  since  it  will 
nj 

always  be  multiplied  by  a  field  which  is  identically  zero  in  a  neighbourhood 
of  the  origin,  the  singularity  can  be  effectively  ignored  in  the  subsequent 


(17) 

(18) 
(19) 


analysis 


,  3uk!*ni„  ,  Uni  2AQ  f  3  ,  3*nj>. 

4  (o)  U  3^  3^  **  *  4  niUuk  3^  °  ^  ‘  4(a) 


3*  3+  ^  3^  . 

*0  "sK  3^f  *  “k  ajp"  “  -  /«.)  \  9^"S^  **■ 


r  3  ,  nl  . 

+  4(a)  IT  {vak)  JZr  ** 

k  i 


after  two  applications  of  the  divergence  theorem  and  Lemma  1.  Since  4*n 
vanishes  when  R  *  o  advantage  may  be  taken  of  (13)  to  assert  that 


l&L  ■  -vK?1] 


*.a--V2R+1)-^ 2 

0 


Inserting  (25)  in  (24)  shows  that  there  is  no  contribution  from  0. 
Hence,  combining  (22)  -  (25),  we  obtain 


iijdx 

r  ~ 


wherein  MY  be  replaced  by  -u>  pu^  by  virtue  of  (2). 

A  further  formula  is  helpful.  It  originates  from 


3Tik*V 


[  IK  n]  r 

JT(o)  3x,  3x  ~  niTi 


nl  _2 


<fdQ  -  L 


*i  3xk  0  1  ik  3xx  T(o)  ik  axi3xfc 


dx  (27) 


after  drawing  on  Lemma  1  again.  Now,  from  (21)  and  the  symmetry  of  the  double 
partial  derivative. 


4 


a2*  au±  a2* 

4(0)  Tik  ^  ”  4(c)  2U  3xk  ax^dx^  ^ 

a2*. 


*2* 

la  2\““i  i^; 02,111  -  4(oi  2ui  ^  **■ 


la  2<V 


a2*  .  .  a*  . 

nj  3u  nj  -)  2  n 

ni3^  a^J"  i°  dQ 


ax,  a 


iaxk 


f  *%  3  r  3u  \  rnj  . 

4(0)  3x,  ^Ui  aXjJ  3x^  ~ 


Also 


3*  .  3u  3u .  du.  a*  . 

4  Vik  ai^  ^  ”  4V  aiT  +  +  a^}  a^  da 


and  in  here  we  can  put 


3uk  3*ni 


ni  ix  3x.  “  %  )x,  9x 


K**} 

i  3xi 


through  (25).  Hence  (27)  -  (29)  supply 

3u 

.'♦2t, 

B  i 


3TiIc  3*ni  ,  3Ub  3  3+n-S  2 

t(o)  a^”  axj^  "  4{x  ax“  *  2  1x^  (wui)}nk  ax^  °  dfl 


(28) 


(29) 


+  24{ni  a^  ‘V  a^f  -  "k  *T±  ■5^i)>°2dQ 


”  4(o2  axi  ^ui  axj^  Jk^  ^  • 


Since  3  ♦jjj/a^3’^  "  0  the  integrand  of  the  second  surface  integral  can  be 
expressed  as 


{  a  a  ^  3*m 

K  a^  *  "k  a^K  a^  * 

But,  if  in  (12)  C  is  allowed  to  contract  down  to  a  point  in  such  a  way 
that  S'  becomes  all  of  S,  we  see  that  the  integral  of  n  a  grad  (  over  a 
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closed  surface  is  sero.  Hence  the  second  integral  over  Q  disappears  and 


2n+1 

; 

o 


3u 


,  4(X  IT  (<’1  +  2  aT  lwl>,snja0 


-  2 1^  (“i  |jH>  5^ 


(30) 


Note  that  in  obtaining  (26)  and  (30)  no  derivative  or  combinations  of 
derivatives  of  other  than  those  occurring  in  Lemma  1  have  been  called 
for. 


A  rewriting  of  (26)  and  (30)  is  desirable.  Put 

Sn)'*'  ’  4  wS<0>Snjdn  - 

anjlo)  ‘  /n  l»4<o)8n:jafl  , 


V'  ‘  4  ‘IT  “."’ls„3dn  ' 


V”  ’ 

W”  *  4‘9 »>s„Jaa 
■  4  ,U1  «r*fl  • 


In  addition  (A.1)  and  (A. 2)  of  Appendix  A  may  be  quoted  to  give  when  n  >  0 


grad  ♦ 


n+1  ,  n-1 

„n+2  _n-1 


nj 


2n+1 


aw-n-1  '  '  ~  &■- 

where  the  summation  convention  does  not  apply  to  a  repeated  Greek  suffix. 
Then,  if  n  >  0,  (26)  furnishes 
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f 


# 


* 


,  . 1  .  n+2  ,  n+1 

■jT  8nj(a)  -  -  ^r)s„j(R>  +  5  i  feoj{Vl,a(R)  +  Vl,«(R>} 

O  R  O  R  <F*-n- 1 


n+1  n-1 

^2n+T  ^n+1  *M{Vl,0(R>  +  *n-1,0(R)} 


(31) 


n+1  n+1  n-1 

“In  i  (-aj)k^n+1,o  k(R)  +  2n+1  ^  (*0j)k^n-1 ,0  k(R)ldR 

R  a— n- 1  o 


while  (30)  provides 


2  n+1 
a" 


n+1 


^  {r  (a)  +  2s  (0»  -  /Jl^  l  CR>  -  2v  (R) ) 

i  R  QF-n- 1 


n+1  n-1 


(32) 


"  *2n+1  ^^,*03  *  ^®nr-1 , 0(RJ  "  2Vl  ,0(R)  }]dR 


for  j  -  -n,  -n+1,...,n.  If  n  ■  0  (which  entails  j  *•  0)  the  terms 
involving  a. .  are  missing  because 

~  p  j 

1  2 


*00  "  \  6odS1,/R  * 

OF—  1 

With  the  convention  that  the  terms  in  a.,  are  to  be  removed  when  n  is 

—  p] 

placed  equal  to  zero,  (31)  and  (32)  can  be  retained  for  all  n  >  0. 

The  notation  will  be  further  simplified  by  putting  (31)  and  (32)  in  the 

form 

n  _  A, (R)  n+1 

*nj(0>  "  2n+1  M  n  +  2n+1  5nj(R)}dR  » 


(33) 


C  , (R)  n+1 


„  C.(R) 

-2—  /*{- SJ -  + 

3**+1  Jn  1  n 


R 


where 


rnj*0^  +  2snj*°*  2n+1  J0  *  Rn  ’  ff2n+1  ‘'nj 


n+1 


D  . (R) }dR 


*nj*R>  "  _R^nj  +  ^~aj*rn+1  ,a  +  *n+1,a*  ~  ^aj^k^n+l  ,ck^ 


(34) 


(35) 


n-1 


Wu  ’ Vi..1  -  'VkVi.ft1 


9 


(36) 


(37) 


T 


n+1 


Cnj*R^  "  ^  ~oj  ^1,a  2~n+1,a*  ' 

J  a—n-1  J 

n-1 


(38) 


The  body  is  of  finite  size  end  therefore  there  must  exist  positive  finite 
o1,o2(o2>o1)  such  that  the  field  under  consideration  is  identically  zero 
inside  the  sphere  of  radius  and  outside  the  sphere  of  radius  o ,  the 

centres  of  both  spheres  being  at  the  origin.  Obviously,  can  be  reduced 


and  o2  increased,  if  convenient,  so  long  as  both  remain  finite  and  positive. 

Since  g^to)  “  0  for  o  >  02  the  right-hand  side  of  (33)  aust  vanish 
for  all  such  a.  However,  the  integrand  also  disappears  for  R  >  and 

accordingly 


V{&nj(R)/Rn}dR  *  42  Rn+\)lR)SR  “  2.  * 


Similarly 


/02{Cnj(R,/Rn}dR  -  0,  /„2  Rn+1Dnj(R)dR  -  0  . 


(39) 


(40) 


Let  p  be  a  positive  integer.  Suppose  firstly  that  p  >  n.  Then,  from 


(33), 


lEnj<°,|2<~ 


2n 


(2IH-1) 
2n 


*»i2> 


,  ,  1  1  „ 

2p-2n-  -  a  -  -2p 

(2ftf„2  {/0R  'V*)!  <» 


2o 


,o  2n*2**  f  1 

+  In* 


dR  f0  R 


a  „2  ”2p  -4n-2,_  .2 


|Bni(R) |  dR) 


by  Schwarz's  inequality.  Hence 
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(41) 


♦  J®  R2  2P|Bnj(R)|2dR}(2rrt-1)"2 


for  p  >  n.  When  1  <  p  <  n  we  draw  benefit  from  (39)  and  write 
A  .  (R)  ,  0  A  (R)  O  2p+  |  -2n  o  -2p-  f 

<  «l  -  l/0  =SV  < /„  *  “/„■>  I4.)'*"  " 

R  R 


for  a  <  o2*  Since 


we  deduce  that 


f°2  2pf  4  “2n  -  2pf  4  "2n 

I a  R  dR  <  ^  R  d 


2P*  a  0,  -2P~  7 


(2n+1  )2  |jgnj  (c)  J  2  <  80  4  /fl2  R  4|Anj(R)|2dR 


_  1  1  . 

2p+  r  _  r  -2p 

+  20  4  R  ISnjtR)!  « 


for  1  <  p  <  n  and  o  <  o2» 

From  (41)  and  (42)  follows 


°»  n  ,  p  n  2p+  —  —  -2p 

l  1  l£_.,<0)|2  <  l  l  4o  fo  R  |A  (R)|2dR(2nt1)' 

»-0  j— n  -1  n-0  j— n  J 

*  *  2pf  4  f°2  "2p“  4  2-2 

+  I  I  80  /  R  |A  (R)|  dR(2n+1) 

‘  n-p+1  J— n  0 

+  I  I  20  2  /?  R2  |B  (R)|2dR(2nt1)’2 

n-0  j— n  0 


for  o  <  Oj. 
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To  estiaate  the  right- hend  aide  of  (43)  observe  froa  (35)  that 


n+1 


I*,!1*"  ,ilJ„>lWV 


.)!' 


n+1 

+  3*  ^  £ttM,ck* 

CF-nrl 


Mow 


n  n+1 

I  I  I  ~aj(rn+1,o  +  1n+1,a)* 
j— n  CF-n-1 


(44) 


.2  n+1 


<n+1K2ifr1.r,  y  |r  +  .  2 

2n+3  ar-n-1  0+1 ' 01  n+1,a|Z 


froa  (A. 8 ) •  Further 
n+1 


3  3  n+1 


I 


£  ,'Vk  V.,*1  ‘ *1  J'J„ , 

■-n-1  J  r-1  a»1  0F-n-1 


and  ao 


n  n+1 


2  3  n+1 


i  I  I 

j— n  CF-n-1  »“1  CF-n-1 


froa  ( 44 ) .  Hence 


n  n  „  _  n+1 

l  l&nj(R)l  "  3°2  ^  l^l2  *  3C2n+-1)2  l  / 

j— n  J  j— n  J  of- n-1 


.  3  n+1  , 

*  «<»■">  I  I  1*^1.  J 


F-1  OF-D-1 


Another  observation  is  that 


I  I  “  /0Ipui2u(R)|2<*« 

n-0  j— n  3 
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1 


« 


«  * 


on  Account  of  the  completeness  of  the  spherical  harmonics.  Since  p  is 
bounded  it  follows  that 


I  X  <  B  /Qia(R)i2dQ  . 

n-0  j— n  J 


Similarly 


•  n+1  _ 

l  l  'Vi,«(R)  *  v,,.(R>'  *  'W  w 

n-0  a— n-1 


J  l  I  'Vi.m'2  <  ■  4'S'*"2" 


•  3  n*1 

X  X  X 

n-0  e-1  a— n-1 

where  B  is  now  being  used  generically  for  a  constant  independent  of  R. 
It  may  be  shown  in  a  similar  fashion  that 


Jjv”'2  *  *S  ♦  3l2n+,,VL,lc-’.R  *  w! 


,  3  n-1  , 

♦  27  (2n+1 )*  l  l  |t  I2 

*-1  0--n+1 


the  last  two  tens  being  absent  when  n-0.  With  this  convention 


—  n-1 


£  „  L'Vi,»  ♦  Vi,#1  *  *  /al,*,R>l  *  'V”.1  >aa 


n-0  0— -n+1 
-  3  n-1 

X 

n-0 


°  n”1  o  1 

K  \  J~, '**-'■»'  <R/a's(R,''“  • 

-0  »-1  0—  n+1 


Let  us  now  suppose  that 

Jq{|jj(R)|2  +  |3uB/3xB|2}dQ  <  CpR2P"  2  .  (45) 
This  statement  is  obviously  valid  when  p  -  1  from  the  assumed  properties  of 
the  field.  Indeed,  the  left-hand  side  is  zero  when  R  <  cr^  and  when 


19 


R  >  ©2«  Then  (43)  and  the  succeeding  Inequalities  reveal  that 


/Q|yu(o) J2dO  -  l  l  |p.(o)|2  <  B  C 
n-0  j»-n  J  p 


2pf  2  2 
a  ‘‘(l+o^) 


for  o  <  o2«  Invoking  (3)  we  infer  that 


2  2  2Pf  2  2 

/Q|u(0)|  dO  <  BK  Cpa  (1+o2)  . 

A  similar  analysis  applied  to  (34)  furnishes 


(46) 


2pf  2 

4,x  ax;  Vff)  +  2  a^  <‘VI2<»q  <BCpa 


(47) 


and  again  (3 ) ,  together  with  (46),  gives 

,  3u  ,2  _  2p+  ^  _ 

4|arl  ‘“V  l,*°2l  • 

Wl 

Inequalities  (46)  and  (48)  may  be  combined  as 

2pf  - 

/  {|u(R)|2  +  I 3u  /3x  | 2}d0  4  BK2C  R  2(1+o2)  (49) 

*•  mm  p  2 

for  R  4  o2» 

Then,  when  (45)  holds  for  some  p,  (49)  shows  that  it  holds  for  p  ♦  1 
2  2 

provided  that  cpfi  *  BK  C^d+Cj) .  since  it  is  certainly  true  for  p  -  1 

Jq{|u(R)|2  +  |3un/3xB|2}dn  4  {BK2(1+02)}P_1C1R2P’  2 
for  R  4  o2  and  p  -  1,2,...  .  An  immediate  consequence  of  letting  p  ♦  • 

1  l 

2  2  2 

is  that  u  =  0  for  R  <  1/B  K(1+o2)  .  It  is  transparent  that  5  0  in 

the  same  region. 
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Let  Rq  “  1/bW2)2.  Choose  another  origin  which  liee  within  the 

sphere  of  radius  It..  The  body  is  entirely  inside  the  original  sphere  of 

e  W 

radius  Og  and  therefore  contained  in  a  sphere  of  radius  +  RQ  about 
x^ .  however,  the  field  is  identically  zero  in  a  neighbourhood  of  x^  and  so 
the  foregoing  theory  aay  be  applied  to  show  that  it  is  identically  sero  inside 
a  sphere  of  radius 

_  1  1 

R1  -  <1  *  (02+Rq)2}  WK  . 

In  view  of  the  arbitrary  selection  of  we  deduce  that  the  field  is 

identically  sero  throughout  a  sphere  of  radius  Rg  +  Rj  centered  on  the 
original  origin.  Clearly,  the  process  can  be  continued  and  the  field  will  be 
identically  zero  inside  a  sphere  of  radius  Rg  +  R^  +  Rj  +•••+  Rjj  ■  Sn  where 

_  i  i 

R  -  {1  +  (S  ,  +  R.)2}  2/R2K  .  (50) 

n  n-l  0 

If  Sn  approached  a  finite  limit  as  n  ♦  »,  it  would  be  necessary  for 

R  ♦  0  which  is  inconsistent  with  (50).  Therefore,  the  body  will  be  totally 
n 

enclosed  by  the  sphere  of  radius  SR  after  a  finite  number  of  steps.  Thus 
the  field  is  identically  zero  throughout  the  body  and  the  theorem  is  proved. 

6.  The  exterior  problem 

The  next  matter  to  be  examined  is  what  happens  when  the  regions  T+ 
and  T_  are  bonded  across  S  and,  instead  of  specifying  the  displacement  and 
traction  on  S,  the  behaviour  at  infinity  is  prescribed.  It  will  be  supposed 
that  the  material  parameters  X,  p,  p  have  in  T+  the  constant  values  , 
|iQ,  PQ  respectively.  The  restrictions  that  will  be  imposed  are  covered  by 
Conditions  A.  In  addition  to  satisfying  (3)  the  material  parameters  must 
comply,  in  both  T_  and  T+,  with  either 

(a)  Xm(p)  >  0,  Im(X+2p)  >  0,  Im(«2p)  <  0  , 


♦ 


< 


# 
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or  (b)  Im(y)  >  0,  Im(A+2y)  >  0,  w2 p  >  0  , 

2 

or  (c)  y>0,  A  +  2|i>0,  «p>0  . 

There  is  no  difficulty  in  checking  that  any  one  of  the  alternatives  is 
sufficient  to  ensure  the  validity  of  the  constraints  placed  on  am  and  afe 
in  (3.  Furthermore,  under  any  of  the  three  conditions, 

lm(iub)i0)  >  0,  Im{iua(A0  +  2 ) }  >  0  .  (51) 

We  shall  now  prove  the  following  theorem. 

Theorem  3.  Let  p  be  continuous,  X  be  continuously  differentiable  and  u 
twice  continuously  differentiable  in  T_  +  S  and  have  the  constant  values 
P0,  Aq,  uq  in  T+.  Let  u,  satisfy  (1)  and  (2)  when  xj^  8  and  be  such 

that  u^  and  3u^/3x^  are  continuous  in  both  T_  +  S  and  T+  +  S.  Suppose 
further  that 

<2>+  -  V+  "  (njTjk)-  (52) 

on  S.  Then,  if  u,  t  satisfy  the  radiation  conditions  (8)  at  infinity 
u  s  0 ,  t  =  0  subject  to  Conditions  A. 

Proof.  Let  R  be  so  large  that  a  sphere  of  radius  R  centered  on  the  origin 
totally  surrounds  S.  Then 


k  njVR)“k(R)R2aD  *  4  nj(Tjkuk)+ds 


p0“  ,a,}d* 


where  T  is  the  volume  between  S  and  the  surface  of  the  sphere  of  radius 

* 

R.  The  star  signifies  a  complex  conjugate.  By  virtue  of  (52),  nj^Tj)cu)c)+ 

* 

may  be  replaced  in  the  integral  over  s  by  n^d^u^) 
divergence  theorem  applied  to  T_.  Hence 


and  then  the 


Iq  njVRlVR)R2d0 


/, 


,3u  ,2  .  3u.  3u.  3u,  3ul 

T(R>U|^r|  +  2  3^*£+  3^ 


)- 


pw2|u|2}dx 


The  Imposition  of  the  radiation  conditions  (8)  now  enables  one  to  say  that,  as 
R  +  «, 


,3u  .2 

. 3u„  3u„ , 2 

.  3u- 

3u_ . 

—a 

—2  4. 

|3x3 

3x2I 

l3u3  3uH2  2  2  ,  2*2 

+  jsx^”  +  3xj|  }  "  p“  l“l  Jdx  -  -/Q{iutow0(|u|  -  lx*u|2)  (53) 

+  iu»a(X0+2p0)  |x«\i|2}R2dO  +  o(1)  . 

Under  conditions  A  the  imaginary  part  of  the  left-hand  side  of  (53) 
cannot  be  negative  whereas ,  on  account  of  ( 51 ) ,  the  imaginary  part  of  the 
right-hand  side  cannot  be  positive.  This  is  possible  only  if 


(54) 


(55) 


,3u  ,2  , 3u.  3u_,2 

. 3u,  3u, . 2 

.  3u. 

3u. . 

Isrl  +“{|3^+3^|  + 

1 — ~  +  — 11 

Ux3  3x2I 

1^ 

»x3l 

-  pw2|u|2]dx  -  0  , 


lim  J0{bu0(|uj*  -  |x*jjl*)  +  a(l0+2p0)|x»u|2)R2da  -  0 


R*« 


It  is  evident  at  once  from  (54)  that,  for  conditions  A(a),  u  =  0.  It 
follows  that  t  S  0  and  the  theorem  is  proved.  Similarly,  under  A(b),  (54) 

implies  that  T  =  0  and  then  (2)  enforces  u  S  0.  Again  the  theorem  has 

been  demonstrated. 

For  conditions  A(c),  (54)  is  no  longer  helpful  and  (55)  must  be  turned 
to.  In  Appendix  B  it  is  shown  that  the  radiation  conditions  entail  (equation 
(B.1 1 ) ) 
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-imaR*  a  e“ia4,R  i 

u  -  -iua  fQ  l  x  -  iutox  a  ^  -  +  0(— ) 

R 


subject  to 

A 

ijg  «x  -  0  .  (56) 

A  A 

Since  x  and  x  a  ^  are  perpendicular  the  implication  of  (56)  ia  that 


|u|2  -  {((ua|f0l)2  +  (atolg^  |  )2}/R2  +  0(1/R3) 


Accordingly, 


Vv^'ao'2  +  (V2“o,u2a3|fo|2}dQ“  0 


stems  from  (55). 

2 

By  hypothesis  the  coefficient  of  Ig^l  is  non-sero  and  so  we  are 

obliged  to  have  gQ  =  0.  Similarly,  fg  =  0.  But,  in  Appendix  B,  it  has  been 

demonstrated  that  ffl  =  0  and  <jg  =  0  make  u  s  0  in  T+.  Hence 

r,  :0  in  T  .  It  follows  that  (u).  -  0  and  (n.T..  )  -  0  on  S.  Prom 

jk  +  ~  +  *’  3  jx  * 

(52),  (u)  -  0  and  (n^t.,  )  -  0  and  now  Theorem  2  enforces  u  =  0  and 

"■  *  "■  j  jk  -  '*  ~ 

T =  0  in  T  .  The  proof  of  the  theorem  is  finished, 
jk  “ 

Evidently,  the  same  method  of  proof  but  drawing  on  Theorem  1  can  be 

applied  to  the  piecewise  homogeneous  body  and  so  we  can  state 

Theorem  4.  Let  T_  be  a  bonded  nested  piecewise  homogeneous  body.  Let  u, 

t  satisfy  (1)  and  (2)  except  on  any  interface  and  have  continuity 

properties  analogous  to  those  of  Theorem  3.  If  u,  comply  with  the 

bonded  boundary  condition  and  satisfy  the  radiation  conditions  at  Infinity 

then,  under  conditions  A,  u  S  0  and  T..  SO. 

- 1 - 1  a  ~  -  jk 

The  theorems  have  been  proved  with  the  boundedness  of  R|u|  as  part  of 
the  radiation  conditions.  However,  it  is  demonstrated  at  the  end  of  Appendix 
B  that  this  requirement  can  be  dropped  without  modifying  the  assertion  that 
f0  s  0,  s  0  implies  that  u  =  0  in  T+.  Therefore,  the  validity  of 
Theorems  3  and  4  is  unaffected  by  this  change  to  the  radiation  conditions.  we 
state  this  as  a  corollary. 
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7.  Uniqueness  for  scattering 

An  easy  consequence  of  the  theorems  of  the  preceding  section  is  that  the 
solution  to  the  scattering  problem  is  unique*  Let  an  incident  displacement 
u*  be  generated  by  some  means*  Zt  is  presumed  that  any  scattered  field 
produced  is  outgoing  at  infinity  and  so  satisfies  the  radiation  conditions* 

If  there  were  two  possible  scattered  fields,  taking  the  difference  between  the 
total  fields  would  eliminate  u*  and  have  a  field  satisfying  our  theorems* 
That  field  must  therefore  be  identically  zero  and  uniqueness  of  the  scattered 
wave  has  been  established. 

Theorem  5.  If  a  given  incident  wave  produces  a  scattered  wave  satisfying  the 
radiation  condition  then*  under  the  conditions  of  either  Theorem  3  or  4,  the 
scattered  wave  la  unique. 

It  is,  of  course,  sufficient  for  uniqueness  to  impose  the  radiation 
condition  in  the  form  (8)  alone  but  it  may  be  more  convenient,  in  practice,  to 
keep  the  boundedness  of  Ru  available  even  though  it  is  jettisonable. 


APPENDIX  A 


This  appendix  is  devoted  to  deriving  soae  properties  of  spherical 
harmonics  which  are  needed  in  the  main  text.  The  summation  convention  is  not 
employed  in  this  appendix. 

Since  Rnsn^(8,+)  is  a  solution  of  Laplace's  equation  so  is 
8(RnSn^ J/dx^.  It  will,  however,  be  of  one  degree  lower  and  hence  is 
expressible  in  terms  of  solid  harmonics  of  order  n  -  1.  Hence  there  are 


constant  vectors  a  ,  such  that 
~pj 


n-1 


grad(R  S  . )  *  ^  I  p  ^  •  (A#t) 

J  p—n+1  'F 

A  similar  argument  reveals  that  there  are  constant  vectors  b  .  such  that 


•  'M1 . .  • 

q»-n  ^ 


(A. 2) 


It  is  understood  that  n  >  1  in  (A.1)  since  SQQ  is  constant  and 

*rad  s00  “  2* 

Prom  (A. 1 ) 


div{R1_2ngrad  <RnS  .>>  "  X  *  i  ’  .  /Rn> 


n-1 


nj 


p^-n+1 


~Pj 


n-t,p' 


n-1  n 

X  •  j  •  I  SqpS  /R1 

p^-n+1  FJ  q— n  q 


rt+1 


on  account  of  (a.2).  Since  V2(Rns  . )  *  0  we  deduce  that 

nJ 

n-1  n 

( 1-2n)nS  «  I  a  .  •  £  b  S 

nj  p«-n+1  q~n  "W 


The  orthonormal  property  of  the  SRj  now  leads  to  the  conclusion 


n-1 

X 

p“-n+1 


»  a  *  *  (1-2n)n5  (q,j  -  -n,...,n)  . 

pj  ~qp  qj 


(A.3) 


26 


A  further  result  of  soee  interest  cen  be  derived 


n-1  n 

P“ 


l  I  ^2n-1)«pi  +  {2“+1)fejp}*{(2n"1)2pj  +  {2n+1)fejp} 

-n+1  j— n  PJ  JF  PJ 


•  f  (2n-1 ) {(2n-1 )n(1+2n)  +  (2n+1 )n( 1-2n) } 
j— n 


+  (2n+1)  l  b  •  «2n-1U  +  (2n+1)b1n> 

p— n+1  j— n  JP 


r «.  J 

~jP 


by  (A. 5)  and  (A.3).  The  first  sun  on  the  right  clearly  vanishes  and  so  does 

the  second  when  (A. 4)  and  (A. 6)  are  invoked.  Therefore 

(1-2n)a  .  ■  (1+2n)b.  (p  “  -n+1,...,n-1»  j  *  -n,...,n)  .  (A. 7) 

~Pj  ~jp 

Zf  (A. 7)  is  incorporated  in  (A.S)  and  (A. 6)  we  obtain 
n-1  .  2 

l 

p*-n+1 
n 

l 

<r 


n-1  2 

slp  ‘  V  '  iia  4jq  <3'IJ  ’  . . nl  ' 

L„  *>q  *  *pq  ~  4>p  . -’>  • 


(A.S) 


(A. 9) 
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APPENDIX  B 


In  this  appendix  soma  properties  of  the  field  in  T+  will  be  obtained 
when  T+  is  a  homogeneous  medium  (with  XQ ,  all  positive  and  finite) 

and  the  radiation  conditions  are  imposed  on  the  field.  Then,  it  has  already 
been  explained  in  §3  that 


"  4{TjiV5'X>ni  “  Vmum  3^  V5'*’ 

-  M0(niuj+njui)  g“  9jk(x.X)}dSy  , 


for  x  in  T+.  Using  the  explicit  formula  for  g^k  in  (4)  we  obtain 


2,2 


u  *  grad  9  +  grad  div  g  +  »b  j 


(B.1) 


where 

19  9  9  -i*l*TCl 

♦'*>  ■  ^7  5^  WirVws  *7  -  y\yvt*  %  • 

i  ,  9  a  «'iub,2nL1 

pjl*>  ■  •  7T7  VVi)-\>VU  ^7  -  vmv  a^1  ■  IjtcI  • 

So  long  as  x  keeps  away  from  8 

V2*  +  «2a2t  -  0  ,  (B.2) 

V2 ^  *  W*b2£  -  0  .  (B.3) 

On  account  of  (B.3),  (B.1)  can  be  rewritten  as 

£  ■  grad  ♦  ♦  curl  curl  £  .  (B.4) 

Mow,  if  |x|  >  |g|,  | x-j^| 1 5l  can  be  expanded  in  a 
uniformly  and  absolutSly  convergent  series  of  powers  of  1/|x|  and  any  number 
of  derivatives  can  be  taken  without  destroying  this  property.  Also  8  is 
bounded  so  that  there  is  a  finite  R'  such  that  8  is  cosqpletely  enclosed  in 

a  sphere  with  centre  at  the  origin  and  of  radius  R' .  With  R,  9,  9  the 

usual  spherical  polar  coordinates,  it  can  now  be  deduced  that,  for  R  >  R' , 


♦(X)  -  e"iWaR  l  f 
^  _  n 

n^O 


(B.5) 


curl  j>(x)  “  e 


-iofcR 


tr*0 


(B.6) 


the  aeries  and  their  derivatives  converging  uniformly  and  absolutely* 

However,  ♦  must  satisfy  (B.2) .  Inserting  (B.5)  in  (B.2)  and  equating  to 
sero  the  coefficients  of  the  various  powers  of  1/R  we  arrive  at 

2i«.(W1)V,  .  (.in  ,  ij!,  i-iis 

sin  0  3  +  (B.7) 


+  n(n+1)f  ■  0 (n-0, 1 ,2, . . . ) 

n 


Dealing  with  each  component  of  curl  g  in  the  same  way  we  have 

3  ^2 

2idlb(nf  1  >V,  .  |j  (.!„  9  -JS 

sin  v  oy 

+  alnt-Vg^  »  0  (n  *■  0,1,2,...)  , 

In  addition,  the  divergence  of  curl  ja  must  vanish.  To  meet  this 


requires 


90  *  i,  “  0  ' 


(B.8 ) 


(B.9 ) 


•iwbSnn  *  *1  +  slrTe  fe  (8in  0  +  k  l*n%>  "  0  ln>0)  (B‘10) 

where  i  ,  i  ,  i  are  unit  vectors  along  the  directions  R,  8,  $  increasing 
respectively . 

By  means  of  (B.4)  -  (B.6)  the  pertinent  expansion  of  u  in  powers  of 
1/R  can  be  written  down  but  explicit  details  will  be  omitted  here.  It  is 


sufficient  to  note  that,  if  f_  5  0,  (B.7)  implies  that  f  =  0  (n>0 ) 

u  n 

because  we  have  assumed  that  ua  ?  0 .  Thus  f  5  0  forces  ♦  5  0  in  R  >  r* , 
Taking  advantage  of  the  analytical  character  of  the  integral  for  we  infer 
that  t  =  0  in  T+.  Similarly,  ^  i  0  implies  that  curl  £  =  0  in  T+. 
Consequently,  we  reach  the  important  conclusion  that,  if  fQ  S  0  and 
3  5  0,  u  5  0  in  T+. 
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On*  useful  formula  Is  that 


-iwaR 


-iubR 


u  -  -iua  fn  ~ 
~  OR 


i.  -iubi.  a  g  - 


*1 - S!  ••  90  r 


+  0(“) 
R 


(B.11) 


REMARK.  In  the  radiation  conditions  stated  in  the  text  and  as  used  above 
it  is  asked  that  Ru  be  bounded  as  R  +  “•  This  extra  requirement  is 
unnecessary  since  (8)  is  sufficient  by  itself  to  justify  all  the  results  that 
have  been  obtained.  The  proof  of  this  assertion  will  now  be  set  out. 

For  convenience  (8)  will  be  repeated  here;  it  is 

R{Vjm  +  i“b‘,0(uj“umxmXj)  +  iwa(  V2VXjVm}  *  0  * 

A 

Multiplication  of  (B. 12)  by  x^  supplies 

RtVjmxj  +  iua(  V2w0)umxm}  *  0  * 


(B.12) 


(B.13) 


The  reality  of  XQ ,  and  pQ  permits  the  observation  that 


4(njTiju!  -  Vijui)R2dQ  -  /s(njTijuI  -  v!jui)ds  • 

In  other  words  the  integral  on  the  left  of  (B.14)  is  independent  of  R 


(B.14) 


although  and  u^  are  functions  of  R.  Now 


i  Vjm +  ^vvvvy +  i“(vji,o)xjV»r 


'Vj-I2  +  (^o,2|YVBxj|2+  {“<V2V}  'V.1  -**bV.(  VVVV 


+  i«{a(X0+2M0)  -  bw0>(uBTjk  -  “.W-fic 


"  'Vjm'2  +  ^o)2^i  -  u.x.xj|2 


{a(X  +2y  )-bp  }  »  2 

*  «(y.2V  '".V]  *  1“l  V2“ol".x.1 


-  (i  -  “v'.i Y2V,|x.Vj'2  *  “2*b“olV2l'ollV.12 


Incorporation  of  this  in  (B.14)  and  application  of  (B.12)/  (B.13)  shows  that 

.2 


4r2(iVjJ2  ‘  'Vj-'j'2  *  ‘,“oi*.V‘j|2,'*(V2V  * 


*  2  *  *  2 

nust  be  finite  as  R  ♦  •».  Since  |x  t.  \d  >  |x  t.  x , |  this  is  possible  only 

m  jB  B  jb  j 

if 

A>|XBTjB<l0xj|2d0'  /Q|XBTjB(R>|2<iQ'  /oluj(R>-u«(R)xmxj|2dQ'  /Q|um{R)XB,2<ia 

are  all  0( 1/R^)  as  R  ♦  •.  Therefore 

/Q|u^(R)|2dfl  <  2  /Q{|Uj-uBXBXj|2  +  lu^x^  |2}d0  -  ,  (B.15) 

R 

^Q|xBTja(R>j2d0  "  0(1/R2)  *  (B. 16) 

The  order  estiaates  (B.15)  and  (B.16)  are  enough  to  enable  us  to  dispense 
with  the  extra  radiation  conditions  since  they  have  been  derived  entirely  on 
the  basis  of  (8)  or  (B.12).  For  the  purpose  of  that  further  condition  was  to 
ensure  that  there  was  no  contribution  froa  the  sphere  at  infinity  in  coning  to 
the  representation  of  u  in  (B.1).  Mow,  the  contribution  of  the  sphere  at 
infinity  to  u(y)  is  given  by 

a  a  , 

•VVijV*'*’  ■  v.“.  V5'l>  ■  VniVVi'  Uf  V;'X>,R  aa 

with  R  >  Ijjf  •  On  the  grounds  of  the  explicit  formula  (4)  g^  can  be 
written 

«jk  “  f JJc  +  hjk 

.[  |i\ 

where  f covers  the  terms  with  e  and  those  with  e  .  But, 

because  of  the  uniform  expansions  in  1/R, 

*  1 

3x  ”  jk  +  0(“2)  ' 

1 

3x  “  "itt4>xihjk  +  0(  2*  ' 

“  R 

as  R  ♦  •  with  |jj|  fixed.  The  integrand  under  consideration  can  now  be 

expressed  as 


# 
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*  A  A  2 

{Vij  +  ^VVVVj  +  iwaMouj}£jkR 

A  A  A  A 

+  {x^T^  +  i«b(  ^o+W0)XaU«Xj  +  iWbWOu j jkR  +  0(lH^  * 

With  the  benefit  of  (9} 

f  jk  *  xjxkFt8,*>e  iWaR/R  +  0(1  /R2 )  , 

*  *  -itifon  2 

hjh  *  ”  ®jJt)H(6#4)e  /R  +  0(1/R  ) 

Hence  the  integrand  becomes 

{xiTiJxj  +  iwa ( \0+2u0 ) xmua}xkFRe"iuaR 

+  {xiTijYk  *  Vik  -  i<JbMo(uk  '  VjXk))HR8’i04>R  +  0(1 Vij1  +  '2l>  * 

The  first  term  need  not  be  discussed  further  on  account  of  (B.13).  As  for  the 
second  it  may  be  rewritten  as 


[{Vijxj  +  lua<V2l,0)umxm}xk  "  Vik  -  ih*WujXjV 

and  its  contribution  therefore  disappears  by  virtue  of  (B.12)  and  (B.13). 
Moreover, 


(/Q|u(R)|dfl)2  <  4*  /a|u(R)|2da-  0(1/R2)  , 

"  0{1/r2) 

on  account  of  (B.15)  and  (B.16).  Hence  the  value  of  the  integral  over  the 
large  sphere  is  zero  in  the  limit  as  R  ♦  •.  Therefore  the  formula  (B.1)  has 
been  verified  solely  under  the  condition  (8).  Since  the  remainder  of  Appendix 
B  is  founded  on  the  representation  (B.1)  alone  its  conclusions  are  still 
correct  without  the  additional  radiation  condition.  Accordingly,  our 
contention  that  (8)  is  a  sufficient  radiation  condition  by  itself  has  been 


vindicated 
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